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' Abstract. We consider a physical system which is coupled indirectly to a 

Markovian resevoir through an oscillator mode. This is the case, for exam- 
ple, in the usual model of an atomic sample in a leaky optical cavity which is 
ubiquitous in quantum optics. In the strong coupling limit the oscillator can 
be eliminated entirely from the model, leaving an effective direct coupling be- 
' tween the system and the resevoir. Here we provide a mathematically rigorous 

treatment of this limit as a weak limit of the time evolution and observablcs 
on a suitably chosen exponential domain in Fock space. The resulting effective 
^ ' model may contain emission and absorption as well as scattering interactions. 

m 
o 

ON i 

O ' 1. Introduction 

vO 

, The motivation for this article stems from the following problem in quantum 

optics, illustrated in Fig.nj Consider the canonical starting point of cavity QED, 
~ Q j. an atomic system in an optical cavity. In many cases such a system is well modelled 

using only a single cavity mode. The cavity is then effectively described by a single 
quantum harmonic oscillator, and the atom-cavity interaction Hamiltonian takes 
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the form 

(1.1) H = E n tfb + E w tf+Eoib + E 00 , 

where Eij are operators acting on the atomic Hilbert space, e\- = Eji, and b, b^ are 
the cavity mode annihilation and creation operators, respectively. Usually one of 
the cavity mirrors is assumed to be perfectly reflective, while the other mirror allows 
some light to leak into the electromagnetic field outside the cavity and vice versa. In 
the Markov approximation P3G], the time evolution of the entire system (consisting 
of the atom, cavity and external field) is described by the unitary solution to the 
Hudson-Parthasarathy 16: quantum stochastic differential equation 

(1.2) dU t = [y/jbdAl - y/jtfdAt -^tfbdt-iHdt}u t , 

where A t , A t are the usual creation and annihilation processes in the external field. 
The transmissivity of the leaky mirror is controlled by the positive constant 7. 

In many situations of practical interest, 7 will be quite large compared to the 
strengths of the atom-cavity interaction. When this is the case, one would 

expect that the presence of the cavity has little qualitative influence on the atomic 
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FIGURE 1. Cartoon illustration of a prototypical problem that falls 
within the framework of this paper. An atom is placed in a single-mode 
optical cavity. One of the cavity mirrors is leaky, thus coupling the 
cavity mode to the external field (which may ultimately be detected). 
In the "bad cavity" limit the mirror is made so transmissive that the 
cavity can be eliminated from the description of the model, leaving an 
effective direct interaction between the atom and the field. 



dynamics: the cavity is then essentially transparent in the frequency range corre- 
sponding to the atomic dynamics, so that the atoms "see" the external field directly. 
Similarly, we expect that measurements obtained from detection of the outgoing 
field (e.g., by homodyne detection) would depend directly on the atomic observables 
and would be essentially independent of the cavity observables. The hope is, then, 
that the time evolution U t can be described in some idealized limit by the unitary 
solution Ut of a new Hudson-Parthasarathy equation which involves only atomic 
operators and the external field, and in which the cavity has been eliminated. The 
goal of this article is to make these ideas precise. 

Previous work. The elimination of a leaky cavity in the bad cavity limit is an ex- 
tremely common procedure in the physics literature — so common, in fact, that most 
papers state the resulting expression without further comment ( "we adiabatically 
eliminate the cavity, giving ..."). Often the equation considered is a Lindblad-type 
master equation for the atom and cavity; in our context, this (deterministic) dif- 
ferential equation for the reduced density operator can be obtained by averaging 
over the field as in [T^]. One method that is used to eliminate the cavity in such 
an equation, see e.g. [2S| r> involves expanding the density operator in matrix ele- 
ments, setting certain time derivatives to zero, then solving algebraically to obtain 
an equation for the atomic matrix elements only. This method is commonly known 
as adiabatic elimination. Though such an approach is not very rigorous, similar 
techniques can sometimes be justified in the context of the classical theory of sin- 
gular perturbations (Tikhonov's theorem [23]). A somewhat different approach, 
see e.g. QT], uses projection operators and Laplace transform techniques. None of 
these techniques are applicable to the question posed here, however, as we wish 
to retain the external field in the limiting model. Hence we are seeking a singular 
perturbation result for quantum stochastic differential equations, which is (to our 
knowledge) not yet available in the literature. 

A naive attempt at adiabatic elimination for quantum stochastic equations is 
made in [7j (see also [S||2Z]). These authors use the following procedure: 

• First, they obtain Heisenberg equations of motion (in ltd form) for the 
cavity annihilator b t = UjbUt and also for the relevant atomic operators. 
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• Next, they set bt = (where the right-hand side is interpreted as "quantum 
white noise" ) and solve algebraically for bt ■ 

• Next, they plug this expression into the atomic equations of motion. 

• Finally, they interpret these equations as "implicit" equations |25| (a for- 
mal analog of Stratonovich equations) and convert to the "explicit form" 
(a formal analog of ltd equations). The latter are considered to be the 
adiabatically eliminated Heisenberg equations of motion for the atomic op- 
erators. 

Attempts at justifying this procedure run into a number of seemingly fatal prob- 
lems. Forgoing the issue of the mathematical well-posedness of "quantum white 
noise" , the approximation 6^ = seems incompatible with the fact that the right- 
hand side is formally infinite. Next we have to deal with the interpretation of the 
resulting equations; even in the classical stochastic case, it is known that adiabat- 
ically eliminated expressions need not be of Stratonovich type (see ^U] for some 
counterexamples); the singular limit is rather delicate and the resulting outcome 
depends on the way in which the limit is taken. Ignoring even this issue, it should 
be pointed out that the implicit-explicit formalism introduced in |25j (essentially 
along the lines of McShane's canonical extension ^lEj) does not even capture 
correctly the ordinary Markov limit in the presence of scattering interactions; com- 
pare the expressions in to the rigorous results obtained in . It is thus highly 
remarkable (if not miraculous!) that we can essentially reproduce the result of 
using the methods developed in this paper (see example 14.41 in section 0J. 1 

Statement of the problem. Following ^Uj, we seek a "method by which fast vari- 
ables may be eliminated from the equations of motion in some well-defined limit. " 
Which limit to take is not entirely obvious at the outset; for example, the naive 
choice 7 — ► oo only yields trivial results (the cavity is forced to its ground state 
and the atomic dynamics vanishes). To define a nontrivial limit, we introduce the 
scaling parameter e > and make the substitution b i— > e~ x / 2 b in (f 1 . 1|) and l|1.2JI . 
The limit e — > + then has the character of a central limit theorem, and provides a 
nontrivial result in which the cavity is eliminated. (Note that similar scaling limits 
are used in projection operator techniques for master equations 

Our approach, then, is to proceed as follows. First we make the above substi- 
tution. Next we switch to the interaction picture with respect to the cavity-field 
interaction. This gives rise to an interaction picture time evolution in which the 
atom is driven by a quantum Ornstein-Uhlenbeck process. The limit e — > + corre- 
sponds essentially to a Markov limit of this equation, and consequently our proofs 
borrow heavily from the methods developed to treat such limits (particularly from 
the estimates developed in |H|). However, our limits are of a somewhat stronger 
character than those considered in P]|21E] as we take weak limits on a fixed domain 



^We also mention 1241 where some results of |7| are reconsidered. The results are only repro- 
duced, however, at the master equation level; in particular, the quantum noise is not retained and 
the implicit-explicit formalism is not used in those sections where results of |7| are considered. 
Though the naive adiabatic elimination procedure used in |3 El 127) is never well-justified, the 
application of the implicit-explicit formalism is particularly suspect in the presence of scattering 
interactions in view of the discrepancy between 1251 and 1141 . Note also that the manipulations in 
|7) rely on the simple commutation relations between position and momentum; they do not work 
at all, e.g., if E a i, are functions of angular momentum operators. It is thus quite surprising that, 
but unclear why, a reasonable answer is obtained in the particular case considered in 0. 
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in the underlying Hilbert space, rather than "limits in matrix elements" where the 
domain depends on e. We also consider, aside from the time evolution unitary and 
the Heisenberg evolution of the atomic observables, the limiting behavior of the 
output field operators (which can be observed e.g. through homodyne detection). 

For concreteness, we will restrict ourselves to the model described by and 
(|1.2(l . This model is already very rich and widely used in the literature in various 
scenarios. Our results can also be extended to more complicated setups, in par- 
ticular to the case of multiple external fields and oscillators along the lines of ; 
the subsequent extension to thermal and squeezed noises is then also straightfor- 
ward through the usual double Fock space construction, see e.g. [T^J, at least in the 
absence of scattering interactions (En — 0). 



2. Preliminaries 

Throughout this article we work on the product Hilbert space f) = f) sys ® f) osc <g) 
f) rcS v consisting of a physical (e.g. atomic) system f) sys , a quantum harmonic oscil- 
lator rj osc = r(C) (describing e.g. a cavity mode), and an external Bosonic resevoir 
f) rcS v = r(L 2 (M + )) (describing e.g. the electromagnetic field). Here r(rj') denotes 
the symmetric (Boson) Fock space over the one- particle Hilbert space f)'. We use the 
following notation for Fock space vectors: |0) S r(f)') denotes the vacuum vector, 
|/) G r(h') denotes the exponential vector corresponding to / S t)', and £ C r(f)') 
denotes the linear space generated by the exponential vectors (the exponential do- 
main). We will also use the subscripts |/) OS c or |/) reS v, and similarly £ OS c) £resvj 
wherever confusion may arise. 

We define the following standard operators: b and ¥ are the creation and anni- 
hilation operators on f) osc , and A t , A\ and A t are the usual annihilation, creation 
and gauge processes on f) resv , respectively ^Hj • We denote the ampliations of these 
operators to f) by the same symbols. For any / £ L 2 (M. + ) and for any real, bounded 
g e L°°(K + ) n L 2 (R + ) we also define the field operators [S] 



A(f) 



f(t)*dA u A(/)t=/ f(t)dA\, A(g)= g(t)dk t . 

Jo Jo 



We recall that the exponential domain £ rosv can be extended to £ I ' CSV D £ ICSV in 
such a way that £^ esv is invariant under the action of A(f), A(f)^ and A(g), see 
e.g. ^3 pp. 61-65]. Similarly £ osc can be extended to £^ sc so that the latter 
is invariant under b, b^ and b^b. This means in particular that the domain T>' = 
f) sys ®£o SC ®£ r ' osv C f) is invariant under finite linear combinations of operators of the 
form E® {b, , tfb}® {A(f), A(fY,A(g)}, and that commutators of such operators 
are well defined on V . Here E is any bounded operator on () sys and (§) denotes the 
algebraic tensor product. We recall also the useful identities |a) OS c = exp(a6 t )|0) O sc 
and |/) rcsv = exp(A(/)t)|0) rosv H3 sec. 11.20]. 

The starting point for our investigation is the rescaled version of 1)1. 1|) and l|1.2fl . 
We consider the unitary solution U t (e), given the initial condition Uq(e) = I, to 
the Hudson-Parthasarathy quantum stochastic differential equation (QSDE) 

(2.1) dUt(s) = [a/J bdA\ - J J &t d A t - T tfhdt -iH(e) dtX U t (e), 
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where the Hamiltonian H(e) is taken to be of the form 

(2.2) H(e) = - Entfb + -= E w tf + -= E Q1 b + E 00 . 

Here Eij, \\Eij\\ < oo are (the ampliations of) given bounded operators on f) S y S , 
Ejj = Eji, and j,e > are positive constants. The existence, uniqueness and 
unitarity of the solution of Q2.1I1 are established in [5] . 

The interaction picture. We are interested in the limit e — > + . As the oscillator- 
resevoir dynamics becomes singular in this limit, the first step we take is to remove 
these dynamics by going over to the interaction representation. To this end, define 
the oscillator-resevoir time evolution VJ (e) as the unitary solution of 

dV t {e) = LjlbdA\- JltfdA t -l 6 + 6dt| V t (e), 

where Vo(e) — I. Existence, uniqueness and unitarity are again guaranteed by UJ. 
We wish to consider the unitary 

(2.3) U t (e) = Vt(e?U t (E). 

Using the quantum Ito rules I16L we find that Ut(s) is given by the solution of the 
Schrodinger equation 

(2.4) = -i f t (e)U t (e), U (e) = I, 
with time-dependent interaction Hamiltonian 

(2.5) f t (e) = V t (e)^H(e)V t (e). 

Note that Vt(s) commutes with any system operator E (on f) S ys), so we have 
E t (e) = U t (e)^EU t (s) = U t (s)^EU t (e). Hence in order to study limits of the 
form lim £ „,o+ E t (e) it is sufficient to consider Ut{s) rather than Ut(e). 

Quantum Ornstein- Uhlenbeck processes. It is convenient to introduce the quan- 
tum Ornstein-Uhlenbeck (O-U) annihilation and creation processes 

(2.6) d t (s) = -M V t (eybV t (e), a t ( £ )t = V t (e)WV t {e). 

V 4e V 4e 

Using the quantum Ito rules, we find that 

T T it 

da t {e) = —dA t - — a t {e) dt, a (e) = - J — b. 
le le V 4£ 

Solving explicitly, we obtain 



= JZ p -lt/2e 



(2.7) g t ( e )=^_, 



J / e lu/2e dA u 



This allows us to express the interaction Hamiltonian T t (e) in the form 

(2.8) f t (e) = - E 11 ~a t {e)^~a t {e) - E 10 h t { £ y - E 01 ~a t {e) + E m . 

7 V7 V7 
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Lemma 2.1. The 0-U processes satisfy (on T>') the commutation relations 

(2.9) [a t (e),a s (e)]=0, [a t (s),a s (s)^] = G s (t - s), 
where the correlation function G e (r) is given by 

Proof. Recall [201 sec. 11.20] the commutation relation on V 

poo 

[A(f),A(g)^}= / f(t)*g(t)dt. 
Jo 

Hence we obtain 

e-^ 2e f e^ 2e dA u ,e-^ 2£ f e^ 2e dA\ t 
Jo Jo 

e _ 7 ( t+s )/ 2£ / e j U /e du _ £ e - 7 (t+s-2tA S )/2e _ £ e _ 7 ( t + s )/2 £ 
Jo 7 7 

where t A s = min(t, s). Now note that t + s — 2< A s = |t — s|. Hence writing out 
the full commutators and using [6, frf] = 1, the result follows. □ 

The function G e (-) has the property of being strictly positive, symmetric and 
integrable with f_ G £ (t) dr = 1. In the limit e — ► + , G s (-) therefore converges 
in the sense of distributions to a delta function at the origin. This means that the 
Ornstein-Uhlenbeck processes at(e), at(e)^ formally converge to quantum white 
noises as e — > + . Let us remark that these processes may now be written as 

(2.10) 5*(e) = 2 / G e {t - s) dA s ~ x [^G E {t) b. 

Jo V 7 

3. Strong coupling limit (e — > 0+) 
For any g £ L 2 (R + ) we define the future and past smoothed functions 



! f g{t + t)G £ (t) dr = 2 f 
Jo Jt 



g+{t,e) = 2 g(t + T )G e (T)dT = 2 g( T )G e (t - t) dr, 



g-(t,e) = 2 f g(t~T)G e (T)dr = 2 f g(T)G £ (t-T)dr. 
Jo Jo 

We will encounter such functions repeatedly in the following. If g were a continuous 

function, we would have the limits 

(3.1) lim g ± {t,e)=g(t). 

e^0+ 

The space L 2 (R+) is much too large, however, to ensure that the e — * + limits of 
the smoothed functions are well behaved; consider for example a square integrable 
function with oscillatory discontinuity (e.g. sin(l/x)). To avoid such unpleasantness 
we will restrict our attention to the set of regulated square integrable functions, 
following |12|. 

Definition 3.1. Let L±(tt+) C L 2 (K+) n L°°(K + ) denote the set of square inte- 
grable bounded functions / on the halfhnc such that the limits /(t^ = lim s ^ 4 ± f(s) 
exist at every point t E R + . We denote by £^ sv C £ ICS v the restricted exponential 
domain generated by exponential vectors with amplitude functions in L±(M.+ ). 
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Before moving on, we make the following remarks: 

(1) Any g G Lj_(R_|_) has at most a countable number of discontinuity points 
(see e.g. chapter 3]). Hence for such g Eq. 1)3. 1[) holds for (Lebesgue-)a.e. 
t e R+. 

(2) Note that if g G L±(R+), then X[o,s]9 € ^±( R +) for am/ s € R+. Hence 
£j^ sv is a suitable choice for the restricted exponential domain used in 
the construction of the Hudson-Parthasarathy stochastic integration theory 

ng. 

For future reference, we collect various e — > + limits in the following lemma. 
Lemma 3.2. For any g 6 Lj_(R + ), the following hold: 5 ± (t,e) > S'(t ± ), 

lim f GM - s)g(s) ds = lim 2 f GJt - s)g + {s.e) ds = 9 ^ + 9 ^ \ 

e^O+Jo ' ^0+ J Q 2 

Moreover, all these expressions are equal to g(t) for (Lebesgue-)a.e. t S R+; hence 
it follows that g^(-,e) ► <?(•) in L 2 (R< ). etc. 

Proof. The hrst statement follows from 



g + (t,e) = j Je-W 2 , (t + £r)dT ^ 9(i+)i 
where we have used dominated convergence to take the limit. Similarly 

^e^ T/2 g(t~ -■r\^ ... . It',, It — 

The third statement follows directly as 



g~(t,e)= I je-^g(t-er) X[ o,tM(r)dT^^g(t-). 



G e (t — s)g(s) ds = - — — - ^" - — — - . 



To prove the next statement, note that 



dr. 



pi pOO pt/\T 

2/ G e (t- s)g + (s,e)ds = / g(r) 4/ G e (t - s)G e {s - r) ds 
Jo Jo I Jo 

But straightforward calculation yields 

4 G e (t-s)G e (s-r)d S = G s (t-T)-G s (T) expf-g) , 

and the result follows directly. Finally, the last statement of the lemma follows 
from the fact that g has at most a countable number of discontinuities, together 
with the dominated convergence theorem. □ 

For / 6 we define the smeared field operator 

/>oo 

(3.2) A(f,e)= / f(t)*~a t {e)dt. 



It is straightforward to obtain the commutation relation 

pOQ -1 pOO 

[A(f,e),A(g,e)^= f{t)*G e g(t)dt = - f(t)*(g+(t,e) + g-(t,e))dt, 
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where we have written G e g(t) = / °° G e {t — s)g(s) ds = ^g + (t, s) + \g (t, e). Using 
(I2.1U[) . we may also express the smeared field as 



(3.3) 



A(f,s)=A(f+) 



f + (0,s)*b, 



where f^~(t) = f + (t,e). The second term ought to be negligible in the e — > + 
limit and indeed, if tp € f) OS c is a vector with \\btp\\ < oo and ii e(g) = \g} T osv is an 
exponential vector (g £ L 2 (R + )), then 



\(A(f,e)-A(f))<p®e(g)\\ \\A(f t 



\\<p®e{g)\\ 



Ml 



\\e(g)\\ V 7 II v 



(f + (t,e)- f(t))*g(t)dt 



V 7 IMI 



More generally, we can define the smeared Weyl operators 
(3.4) W(f, e) = exp e )t - A(f, e)} , 

which satisfy the smeared canonical commutation relations 



(3.5) 



W(f,e)W(g,e) = W(f + g,£) exp 



-Urn I f(t)*G E g(t)dt 

10 



The smeared Weyl operator W(f, e) ought to converge to the standard Weyl unitary 
W(f) = I osc ® exp{A(/) t - A(f)} ase^ + . It is indeed not difficult to establish 
that for arbitrary /i, • • • /„ € L±(R+) 



||(W(/i,e) •••#(/«,£) 



•W(/ n ))^®e(ff)|| -^1*0. 



This time, no restriction needs to be placed on tp. This is a type of quantum central 
limit theorem, however, it is less abstract than the "limit in matrix elements" 
traditionally encountered in the quantum probability literature ^ |21 114| since the 
limit is taken on the fixed domain f)osc®£rcsv The limiting operator is thus defined 
on the same Hilbert space, though it acts non-trivially on the noise space only. 



4. Limit dynamics 

The limit of the process {Ut(s) ■ t > 0} as e — > + is reminiscent of the Markov 
limits that have been widely studied in mathematical physics 0121 El- Comparison 
with previous results suggests that the limit be again described by a quantum 
stochastic process {U t : t > 0}. We wish to deduce this limiting process by studying 
the limit of matrix elements (tpi\Ut(s)tp2) for arbitrary vectors of the form 



ip i = v i ® \ai) osc <g) \fi] 



Vi G t) sys , a t e C, fi e L\{ 



In other words, we would like to obtain U t as the weak limit of Ut(s), as e — > + , 



on the domain l) sys <8>£, 



syswosc^rcsv ■ 



We will similarly study weak limits of observables 



Ut(e)^ EUt(e), where E is a system observable, on the same domain. 

Formally, we may expand Ut{s) as a Dyson series (by Picard iteration): 



(4.1) 



U t (e)=I + Y(-i) 



71=1 



d.Sr. 



A„(t) 



■dsif Sn (e)---f Sl (e), 
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where the multi-time integrals are taken over the simplex 

A n (t) = {{s n , . . . , si) : t > s n > ■ ■ ■ > si > 0}. 
The Dyson series expansion of (tpilUtis)^} is given by 

(4.2) (Vi|^ 2 ) + Yl-iT / ds n ---d Sl (ij 1 \f Sn (e)---f Sl (e)^ 2 ). 
n =i ■ / A„(t) 

The usual existence proof for differential equations by Picard iteration suggests 
that the Dyson series are convergent. Following m21E], our basic approach will 
be as follows. First, we obtain an estimate of the form 



ds n ■ ■ ■ dsi (^i|T Sn (e) • • • T SI (e)i/} 2 } 

A„(t) 



where is independent of £ and such that X^^Li ^( n ) < °°- This establishes 

uniform convergence of the Dyson series l|4.2|l on e > (by the Weierstrass M-test). 
Consequently, we may exchange the limit and the summation in the Dyson series: 
i.e., we have established that (ipi\U t (e)ip 2 } converges as e — > + to 

oo „ 

(4>M + Y,(-i) n lim / ds n ---d Sl (^ 1 \f Sri (e)---f Sl (e)^ 2 ). 

It then remains to determine the limiting form of every term in the Dyson series 
individually. Summing these we obtain an (absolutely convergent) series expansion 
for the limiting matrix element, which we identify as the Dyson series expansion 
of the solution Ut of a particular quantum stochastic differential equation. This 
completes the proof. Details can be found in section |SJ 

In principle we should establish the results sketched above for every pair of 
vectors ipi. It is convenient, however, to reduce the problem to the study of the 
vacuum matrix element only. To this end we use the following identity: 

ipi = exp-ja,^} exp{A(/ i ) t } v, <g> |0) osc ® |0) reS y 

By commuting the operators exp{ai&^} and exp{A(fi)' f } past T s „(e) • • • T Sl (e) we 
can express the matrix element in the integrand in terms of the vacuum, provided 
we make some simple modifications to T t (s). Hence we have to go through the 
proofs only once using the vectors V{ <g> |0) osc ® |0) reS v 

Lemma 4.1. Define $ = |0) osc ® |0) resv . Then 

{^x\U t {e)^ 2 ) = (vi®®\U t (E)v 2 ®$) (ai|a 2 ) OS c (/i|/2>resv, 

where Ut(s) is the modification of Ut(s) obtained by the replacements 

/4e 

5 t (e) i — ► a t (e) = a t (e) + f 2 (t,e) - W — G e {t)a 2 , 

V 7 

V 7 

Proof. Key here are the simple identities (on V) 

be ^ = e ab ' (b + a), A(g)e A ^ = e A ^ (a( 9 ) + J™ g(s)* f(s) ds^j , 
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which allow us to write using p.lOjl 

a t (e)e ab 'e A ^ = e^V^ (a t (e) + f~ (t, e ) - ^ G e (t) a) . 
Hence starting from any term in the Dyson series of the form 
(•0i|T s „(e) • • • f Sl ( £ )V> 2 ) = (vi ® $|e^ b e A(/l) T fl „(e) • ■ ■ f Sl ( £ )e A(/2)t e Q2fc+ i* ® 
the result follows using the above relations if we use additionally that 

e ^ e A(/ l)e A(/ 2 )V 26 t = e A(f^ e ^ ea lb e A {h ) (ai | a2 ) osc (/l |/ 2 ) rcsv 

and that e < b e A ^ $ = $. □ 

The Dyson expansion for the matrix element {ipi\Ut{e)ip 2 ) may now be written, 
up to a constant prefactor of (ai|a2) sc (/i|/2)resvj as 

oo „ 

(vi\v 2 ) + V(-!)" / ds w -.-dai<i;i®*|f. w (e)..-f. 1 (e)«2®*). 

„=1 JA n (t 



Here T t (e) is obtained from T t (e) by making the translations above, that is, 

(4.3) T t (e) = -E ll at{e)a^{e)-^-E lQ ai{e)-^-E 01 a^{£) + E 00 . 
7 V7 n/7 

As the new processes (e) are linear in the original O-U processes, we may write 

T t (e) = E^{t,e)at{e^a t {e) + E 10 {t,e)a t {e)^ + E 01 (t,e)a t {e) + E 00 (t,e). 

The coefficients Eij(t,e) are easily worked out, however, our main interest will be 
in their limit values: we have for (Lebesgue-)a.e. t £ R+ 

£ n (t, £ ) =Mki(t) = -E n , 

7 



lim 


E 10 {t,e) 


= E w (t) 


2 


Eio ■ 


+ -E n f 2 (t), 

7 


lim 


E i(t,e) 


= E 01 (t) 


2 


Eqi 


+ -h(tyE n , 

7 


lim 


Eoo{t,e) 


= E 00 (t) 


= E 00 — 


2 


E m f 2 {t) - 4= 

V7 



0*^io + -/iW*£ii/ 2 (t), 

7 

the limits being uniform in the strong topology. Note that these limits depend only 
on the functions fi describing the resevoir: the parameters a, for the oscillator have 
disappeared, indicating that the oscillator is indeed eliminated as e — > + . 
The expansion of (ipi\U t (e)if) 2 ) /{ai\a 2 ) oac (/i|/2)resv may now be written as 

oo „ 

(vi\v 2 ) + y^( - *)™ / dsn-'-dsi 

• • • (vi\E an p n (s n ,s) ■ ■ ■ E ai 1 (s 1 ,e)v 2 ) 

x m\as n i^] an [a Sn {e)f- ■ ■ ■ [a si (eft* K (^)f 1 , 

where (3k are summed over the values 0, 1 and we write [x]° = 1, [a;] 1 = x. It is 
this form of the expansion that will be most useful in the proofs (section |HJ. 
We can now state the main results. 



n=l 

x 
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Theorem 4.2. Suppose the system operators Eij are bounded with ||£ai|| < 7/2. 
Then there exists a unitary quantum stochastic process {Ut : t > 0} such that 

lim (V>i|t/tO)^2> = (i/>i\U t ih) 

for any pair of vectors 1^1,2 & f)sys®^osc®^re SV • The process {Ut '■ t > 0} satisfies a 
QSDE of Hudson-Parthasarathy type 

dU t = \ {W -I)dk t + LdA\ -DwdAt ~ -ULdt - iHdtl U t , U = I, 

where the coefficients are given by the expressions 

■ L= T^TTf Ew > H = E 00 + lm{Eoih/2 + tE 11 )- 1 E 10 }. 
j/2 + iEii j/2 + iEii 



Theorem 4.3. Under the conditions of Thm. \4-S\ we have convergence of the 
Heisenberg evolution: for any bounded operator E on t) sys and ipi € f) S ys®£osc®£re S v 

lim {ih\Ut(e)iEUt(e)ih) = lim {ipi\U t {e)^ EU t {e)^ 2 ) = (falty EU t ip 2 ) ■ 

Let us demonstrate these results for some models used in the physics literature. 
Example 4.4. Doherty et al. [7] consider the following system, in our notation: 

7 = 2k, H = E 00 - ^ cos 2 (k L x)b' t b, 

where x is the atomic position operator on fj sys = L 2 (R) and Eqq is a free Hamil- 
tonian. 2 According to Theorem 14.21 the limiting time evolution is given by 

dUt = {<W-I) dA t ~ iE 00 dt\ U t , W = K + l i C0S \l LX) '^ 
L J k — igQ cos {Klx) I A 

provided that ||-E7ii|| = g 2 /A < k. According to Theorem 14. HI and the quantum Ito 
rules, the limiting Heisenberg evolution of an atomic operator X is given by 

d{u}xU t ) = U} (i[E 0Q , X] dt + (W^XW - X) dA t ) U t . 

Compare this expression to Eq. (2.16ab) in 7 , taking into account the identity 
exp{2i tan _1 (a;)} = (cos{tan _1 (a;)} +i sin{tan _1 (a;)}) 2 = (1 +ix)/(l — ix). 

Example 4.5. The following interaction Hamiltonian is often used to describe the 
coupling between a collection of atomic spins (total spin J, i.e. f) sys = C 2J+1 ) and 
a far detuned driven cavity mode (see e.g. |22p: 

H = X F z tfb + U{tf + b)+E 00 . 

Here Eqq is a free atomic Hamiltonian and xM are rea l constants. By Theorem 
14.21 the operators W , L, H in the limiting QSDE become 



j/2 + ixF z ' 7 /2 + i X F,' 7 2 /4 + x 2 i 712 ' 

provided ||J5n|| = xJ < 7/2. A common assumption in the literature (a reasonable 
one if the adiabatic approximation is good) is that ||2x-F z /7|| = 2% J/7 <C 1; the 
conventional adiabatically eliminated master equation (such as the one used in 



2 Technically their Hamiltonian _Eqo = Px/^ m i s unbounded, but we sweep this under the rug 
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|22j ) is now recovered by calculating the master equation corresponding to Ut, then 
expanding L and H to first order with respect to 2xF z /j. 

5. Limit output fields 

Aside from the limit dynamics of the system observables, as in Theorem 14.31 
we are also interested in the limiting behavior of the resevoir observables after in- 
teraction with the system and oscillator. In optical systems, for example, these 
observables can be detected (using, e.g., homodyne detection [3]) and the observed 
photocurrent can be used for statistical inference of the unmeasured system observ- 
ables (quantum filtering theory @JE])- The behavior of these observables in the 
singular limit is thus of significant interest for the modelling of quantum measure- 
ments. 

To investigate the limit of the field observables we study the convergence of ma- 
trix elements of the form (ipi\U t (e)^W(g t ])Ut(e) ip 2 ), where V>i,2 € f)sys®£osc®£ r tsv> 
W(g) is the usual Weyl operator with g G and gt]{s) — g(s)X[o,t]( s )- Note 

that unlike in the system operator case, Vt(s) does not commute with W(g t ])- How- 
ever, we obtain using the quantum Ito rules 

V t (e)^W(g t] )V t (e) =exp{A(g t] y - A(g t] ) -2(A(g t] ,e)^ - A(g t] ,e))} , 

or, expressing this in terms of the usual Weyl operators, 

/4e 

7 

Hence we can write 

(5.1) U t (e)^W{g t] )U t (e) = U t (e?W(g t] - 2g+ e )B( y/teftg+^O, e))U t (e), 

where denote by B{a) = exp{a&^ — a*b} the Weyl operator for the oscillator. Using 
the Dyson series for C/ t (e), we now expand (if>i\Ut(e)^W(gt])Ut(e)i/j2) as 



V t {efW{g t] )V t {e) = W(g t] ~ 2g+ e ) exp [ M ( S J(0,e)&t - 5 +(0, £ )*&) 



OO OO „ „ 

V V H)"- m / dt m ■ ■ ■ dh / ds n ■ ■ ■ d Sl 

n, n JA,„(t) JA,At) 



„=0m=0 JA m (t) JA„(t) 

x(^i|f tl (e) • • • t tm (e)W(g t] - 2g+ e )B(y/4^g+(0,e))f Sn (e) ■ ■ ■ T ai (e)^) 

(for notational simplicity we have used the convention f/± f t \ ■ ■ ■ = I here). The 
limit of this expression is most easily studied by commuting the Weyl operators 
through the T s (e) terms, in the spirit of lemma |4~T1 In particular, using 

W(f) = exp{A(/)t} exp{-A(/)} exp / -i [°° \f(t)\ 2 dt 



and similarly 

B(a) = exp{a6 t } exp{— a*b} exp{— |a| 2 /2}, 

then moving the conjugated terms to the left and the remaining terms to the right 
(where they operate trivially on the vacuum), the problem can be reduced to the 
manipulations used in the proof of Theorem 14. 31 Details can be found in section 
We can already guess at this point, however, what the answer should be. As 
g t ](s) — 2gj(s,e) — > — <?t](s) s-a.e., and as y/Ae/^y g^(0, e) —> 0, we expect that 

U t (eyW(g t] - 2g+jB(y/tefig+(p,e))U t (e) u}w{-g t] )U t . 
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This is in fact the case. 

Theorem 5.1. Under the conditions of Theorem \4-^ we have the following: for 
any Vi,2 € f) sys ®£ osc <8>£ r t, v and g € L±(R+) 

E lim fa\U t (ejtW{g t ])Ut{e)ih) = (il>i\U}W(-g t] )Utih)- 
6. Proofs 

In the previous sections we have set up the problems to be solved, and we have 
investigated in detail the Ornstein-Uhlenbeck noises and the associated correlation 
functions and limits. With this preliminary spade work at hand, the remaining 
(technical) part of the proofs, as outlined in section 0] follows to a large extent 
from the proofs and estimates in |14) . Below we work through the required steps 
in the proofs, however we refer to f° r some detailed calculations. 

6.1. Proof of Theorem IP1 

Wick ordering. All steps of the proofs require us to evaluate the matrix elements 

<$| [a Sn ( £ )t]«- [a Sn ( £ )f " • • • [a st ( e )t]<* [a Sl ( e )]* *) 

that appear in the Dyson series. The solution to this problem is well known and 
proceeds by applying Wick's lemma For given sequences a = (aj), /3 = 
define the sets P(a) — {i : a.{ = 1} and Q((3) — {i : j3i ~ 1}. Let 3(a,/3) be the 
set of all maps J : P(a) — > Q(f3) that are bijections and that are increasing, i.e. 
J(i) > i. Then by Wick's lemma we obtain 

<$|[< i ( £ )t]""[<( £ )]' 5 "...[a Sl (e)t]^[«, 1 (e)] ft <i>)= £ J] G £ ( Sj(l) - 

JeZ(a,f3)itP(a) 

Diagrammatically, this can be represented as follows. Write n vertices on a line: 



S„ Sj s 2 Si 

For each vertex j, draw ingoing and outcoming lines corresponding to the values of 
aj and /3j, as follows: 

aj = ft = 1 a 3 = 1, /3j =0 a ■ = 0, ft = 1 a, ■ = ft = 

:v: .\ r: 

Sj Sj Sj Sj 

Note that P(a) is the set of vertices that have outgoing lines, whereas Q(j3) is the 
set of vertices that have incoming lines. Next, connect every outgoing line to one of 
the incoming lines at a later time (i.e. form pair contractions), in such a way that 
all the lines are connected to exactly one other line. For example: 




S9 Ss S7 S6 S 5 S4 S3 S2 Si 

The ways in which this can be done are in one-to-one correspondence with the 
elements of 3(a, /?); the contracted vertices are then simply the pairs (i, J(i j) where 
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i € P (a). Wick's lemma tells us that the sum over all such (Goldstone) diagrams 
gives precisely the vacuum matrix element we are seeking. 

Step 1: a uniform estimate. Our first goal is to find a uniform (in e) estimate 
on every term M s {n) in the Dyson series: 



M e (n) = 



/ ds n ■ ■ -dsi ^2 ' ' ' E ( v i\Eo< n fS n {s n ,e) ■ ■ ■ E aij3l (si, e) v 2 ) 
x <*|[5. B (e) t ] Q "[5 aB ( e )]' J » ■ ■ ■ {asAe^rK&f 1 *> 



Note that in this expression there is a summation over a and (3. Hence every n- 
vertex Goldstone diagram is going to appear in the sum when we apply Wick's 
lemma, not just those with fixed incoming/outcoming lines for each vertex (as for 
fixed a, (3). Whenever this is the case it is convenient, rather than first summing 
over Z{a, (3) and then over a, (3, to arrange the sum in a slightly different way. 

Every n-vertex Goldstone diagram can be described completely by specifying a 
partition of the set {1, . . . , n}; each part of the partition corresponds to a group of 
vertices that are connected. For example, the nine-vertex example diagram above 
corresponds to the partition {{1, 3}, {2}, {4, 6, 8}, {5}, {7, 9}}. The corresponding 
values of a and (3 are easily reconstructed: a singleton vertex has a = f3 = 0, and 
for a doubleton or higher the first vertex has a = 1, (3 = 0, the last vertex has 
a = 0, (3 = 1, and the vertices in the middle have a — (3 — 1. The sum over 
a, (3 and 3(a,/3), which appears in the expression for M s (n) after applying Wick's 
lemma, can now be replaced by the sum over all partitions <B„ of the n-point set. 

We need to refine the summation a little further. To every partition we associate 
a sequence n = (rij)j^ of integers, where nj counts the number of j-tuples that 
make up the partition (e.g., the example diagram above has n = (2, 2, 1, 0, 0,0,.. .)). 
We denote by E(n) = . jnj the number of vertices in the partition (so E(n) = n 
for a partition of n vertices) and by iV(n) = ^ • rij the number of parts that make 
up the partition. Of course there are many partitions that have the same occupation 
sequence n; the set of all such partitions is denoted 2J n C 23_E( n )- Summing over 
25„ now corresponds to summing first over *B n , then over all n with E(n) = n. 

We are now ready to bound M e (n). First, note that 

\(vi\E anl3n (s n ,s) ■ ■ ■ E aif 3 1 {si,e)v 2 }\ < \\vi\\ \\v 2 \\ C Qrl(3 „ • • • C Ql ft, 

where C a p are finite positive constants that depend only on ||.E a| g||, max[ f ] / 12 
and 7. In particular, C\\ = (4/7)|J_E 11 || and we will write C = max a ^ C a p. For any 
a, (3 corresponding to the occupation sequence n, the number of times that a = 1, 
(3 = 1 will be £ j>2 (j - 2)rij = E(n) - N(n) + n x . Hence 

r< n s- — JV(n)+ni / ^AT(n)-rii 

°a„/3„ ' ' ' Wi/3i S O n O 

We can thus estimate M e (n)/ \\vi\\ \\v2W by 



E(n)=n 

E 



c E W -N^ lcN(a) - ni f ds n ---d Sl TT G s ( Si -sj), 

Jam , ^ „ 



peSn A " w (»,J')~P 



SINGULAR PERTURBATIONS OF QSDE 



15 



where (i,j) ~ p denotes that the vertices i and j are contracted in the partition p. 
A clever argument due to Pule can now be extended to show that 



1 t N ^ 

J e [Si — Sj) < — - 



e / d Sn ...d S1 n G e ( Si - Sj )< n ^ 2L . ini _^ nr 



Essentially, the trick is to rewrite the sum over Q3 n of integrals over the simplex as 
a single integral over a union of simplices, which can then be estimated; see |14l 
section 7] for details. If Cn > 0, we obtain the following estimate uniformly in e: 



Af e (n) < n(n) = |M| \\v 2 \ 



E{n)-n e AE(n)+BN(n) 



^ ni!n 2 !--- 



where A = log(C n /2) and B = \og{t V 1) + log(C 2 V 1) + log(C n 2 V 1) + log 2. 
Summing over n, we obtain 

]_ e AE(n)+BN(n) 00 00 e (kA+B)n f & A+B 



AEla)+BN!n) 00 00 JkA+B)n ( 

n n 1 * I 1 n=C\ *• 



Ml II ^2 II 

' " " " n n " " fc=l n=0 

provided that e A — Cn/2 < 1, i.e. the sum converges provided that ||£ai|| < 7/2. 
Recall that this was a condition of Theorem 14.21 If Cn = we obtain a slightly 
different estimate, which is however even simpler to sum (most terms vanish). 

Now that we have a uniform estimate, the Weierstrass M-test guarantees that 
the Dyson series converges uniformly in e. Consequently, we can calculate the limit 
of the Dyson series as e — > + simply by calculating the limit of each diagram 
independently, then summing all these terms. This is what we will do below. 

Step 2: principal terms in the Dyson series. The contribution of a single Gold- 
stone diagram to the Dyson series has the form 

/ ds n ---ds 1 (v 1 \E anl3n (s n ,e)---E ail 3 1 (s 1 ,e)v2) Y[ G t( s J(i)~ s i) 

for some J € 3(a, /3). A diagram will be called time-consecutive if J(i) = i + 1 
for every i £ P(ct). We claim that in the limit e — > + any diagram that is not 
time-consecutive vanishes: hence we only need to retain time-consecutive diagrams. 
To see this, first note that the magnitude of the diagram above is bounded by 



C n \\vA 



"2 1| / ds n ---dSi JJ_ G e(sj(i) ~ Si). 

-'A„(t) x 



The limit of the latter integral is not difficult to evaluate explicitly. In particular, 
if J is not time-consecutive then the integral vanishes in the limit e — ► + . For 
example, suppose that J(i) =/= i + 1, so that sj(i) > s»+i a.e. in A n (t). Then 

dsi G £ (s J(l ) - Si) ► for any s J(l) > s t+1 

by dominated convergence, as G e (sju\ — Si) is uniformly bounded on [0, Sj+i] when- 
ever Sj(i) > Sj+i and G s (sju\ — st) — > pointwise. On the other hand, 

rn+i rs J{i) j 

/ ds, G e {s J{i) - s^ < / ds l G e {s J{i) - s») = - for any s J{i) > s i+1 . 
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Hence we have by dominated convergence 



I 



dsj(i)---l dsi G e (sj(i) - > 0. 



o Jo 

Proceeding in the same way, we can show that any diagram that is not time- 
consecutive vanishes as e — > + |14l lemma 6.1]. 

It remains to consider the time-consecutive diagrams, for example: 

..:.x:..:. cxzx 

Sg S S S 7 S 6 S 5 S 4 S 3 S 2 Si 

These diagrams have a particularly simple structure: any such diagram is uniquely 
described by listing, in increasing time order, the number of vertices in each con- 
nected component. For example, the diagram above is described by the sequence 
(3, 1, 2, 3). In this way, any n-vertex diagram with m connected components is de- 
scribed by a set of integers ri, . . . , r m such that r± + • • • + r m = n. Now suppose 
that J £ P) is a time-consecutive diagram that is described by the sequence 
ri, . . . , r m with r\ + • • ■ + r m = n. It is not difficult to verify that 



/ ds n ---ds 1 E(si,...,s n ) JJ G e (sj( l ) - s») - 
1 



dt m ■ ■ ■ dt\ S( fl, tx) • ■■;tl , t2 , ti , . . . , ^2 i • ■ ■ i ^mi W ■ ■ : Un ) 

n times r 2 times r m times 

for any function S £ L±(R™ ). Note that n — m is precisely the number of contrac- 
tions in the diagram n, . . . ,r m . 

Step 3: resumming the Dyson series. We now compose the various steps made 
thus far. Starting from the nth term in the Dyson expansion, using Wick's lemma, 
retaining only the time-consecutive terms, and taking the limit as e — > + gives 

/ ds n ---ds\ ^2 ■ ■ ■ ^ (vi\E an n (s n ,e) ■ ■ ■ E aif 3 1 (s 1 ,e)v 2 ) 

• 7a ™W a„/3„ ax/3! 

x (*| [a. B (^P [a. B ■ ■ • [5., (^P [a Sl ( e )]* $) 

T\-\ \-r m —n 



e-*<T 



1 it ~ ' m — '* -i p 

E E ^••■dii(«i|^ ) (i™)---^ (ri) (ii)« 2 ), 

m n,...,r m >l JA m {t) 



where we have written 
Let us now sum all the terms in the limiting Dyson series: this gives 



E 00 {t) r = 1, 

E 01 (t)(E n (t)) r - 2 E 10 (t) r>2. 



r m —n 



k.^)+Eh"E E™ 

n— 1 m n ,. . ,,r m > 1 



x / dt m ---dt 1 (v 1 \E^\t m )---E^\t 1 )v 2 ). 

JA m {t) 
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Now use the fact that n — m = Y^,k( rk — 1) to rewrite this expression as 



/A m (t) » * m ^ " 



v r m >l 



E -7T2^T ) «»>• 



But note that we can sum 



£< r >(t) - , x 1. ., L 

E^T--^ooW-^oiW(E 



">i 



(£n(*)) r 

(2i)' 



Bio(«) 



-^oo(t)-^oift) 1+ ,^ w/2 



^10 (*) 



provided that ||z£Jn(f)/2|| = (2/7)||-Eu|| < 1, which was already required for uni- 
form convergence of the Dyson series. Finally we define 



J a0 



-i E a p — E a i 



1 



7/2 + iE l 



E 



V7 



and note that we can write 

-^-^M l + ^(t)/2 



a/3 



Hence the Dyson expansion for {ipi\Ut(e)ip2} / {ai\a2} osc (/i|/2)resv may be written, 
in the limit e — *■ + , as 



x / dt m • • -dh [h(t m )T m {f2(t m )f™ ■ ■ • [A(*i)T 1 [/2(ii)] /31 . 



Step 4- the limit unitary. It remains to investigate the relation of the limit- 
ing Dyson series given above to the unitary evolution Ut- Consider a Hudson- 
Parthasarathy equation of the form 

dU t = {in dA t + Lio dA\ + L i dA t + L 00 rfij U t . 
By Picard iteration, the solution Ut can be developed into its chaos expansion 

^ J +EE-E/ L -^ m ■ ■ ■ L axPl dkz: 0m ■ ■ ■ dK lPl , 

m a m (3 m oift "' A ™( t ) 

where we have used the Evans notation A^ 1 = A t , A^ = A\, A^ 1 = A t , A^ 10 = t (see 
e.g. ^1 page 151]). Using the usual formula for the matrix elements of stochastic 
integrals, it is evident that (ipi \Ut ip 2 ) coincides with the limiting Dyson series 
above. It remains to notice, as is verified through straightforward manipulations, 
that L n = W - /, Lio = L, L 01 = -DW 7 and L 00 = -iH - DL/2. The proof of 
Theorem 14.21 is complete. 
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6.2. Proof of Theorem 14.31 Conceptually, little changes when we are interested 
in the Heisenberg evolution. Using the Dyson series for Ut(e), we now expand 

(^lUtieyxUtie)^) as 



E E H)"~ m / dtm ■ ■ ■ dh / ds n ■ ■ ■ d Sl 

n n J A™ (t) JA„(t) 



„=0m=0 JA m (t) JA n (t) 

x (^i|f tl (£)••• T tm (e)Xf Sn (£)••• f ai (e)V 2 ) 
(for notational simplicity, we write from this point on J^oCt) *" = I). This equals 

OO OO p p 

n=0m=0 JAm(t) JA n (t) ft m i>„, Ml"l a„^„ Qift 

x (-0i|£' JulWl (ti, e) • ■ ■ E^ mVm (t m ,e)XE an p n {s n ,£) ■ ■ ■ E a ^{sx,e)i>2) 

x [a s „ (e)t]<- [a.„ (e)]^ ■ • • [a Sl ( E )t]* [a Sl $ ) , 

where we have applied lemma ETT1 As before, we can use Wick's lemma to evaluate 
the vacuum matrix element. Drawing vertices on a line in the correct order, assign- 
ing incoming and outgoing lines according to a, /3, /x, v, and connecting them up, 
allows us to represent the vacuum matrix element as a sum over the usual diagrams. 
For example, a possible diagram in this case might be: 



ti t 2 t 3 t& t 5 te t 7 t s tg s 7 s 6 s 5 s 4 s 3 s 2 si 

Note that we do not need to worry about the time ordering (which is obviously not 
satisfied in this case), as the commutators between d s (e) and a t (e)t are symmetric 
in s, t; hence only the order in which the a's and a^'s occur will matter, and we can 
expand in terms of pair contractions in the usual way. 

The first question that needs to be resolved is whether we still have uniform 
control on the convergence of the Dyson series. This does turn out to be the case. 
The argument used previously to obtain the required estimates can be generalized 
also to the Heisenberg evolution, though the details of the argument are somewhat 
more involved in this case. We refer to |141 section 9] for further details. 

The next problem is to determine which diagrams survive in the e — * + limit. 
It is not difficult to see that diagrams with contractions between s-variables which 
are not time-consecutive or between t-variables which are not time-consecutive will 
vanish in the limit; this follows directly from the previous arguments. Hence all 
surviving diagrams must have only time-consecutive contractions within the s- and 
t-blocks. On the other hand, note that we are not integrating over the simplex 
A m +„(t), but rather over the product of simplices A m (t) x A n (t). Therefore con- 
tractions between s- and t- variables do not necessarily give vanishing contributions, 
provided that the corresponding lines in the diagram do not cross — in the latter 
case the contraction would force Sj = ti and Sj = in the limit e — > + , whereas 
integration over A m (t) x A n (t) requires < Sj and tf. < t\. For example, 
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tl tl tz t& t$ t$ t-j tg tg 

must necessarily vanish, whereas the diagram 



S7 Sq S 5 Si S3 S 2 Si 



tl h t 3 tl t 5 t 6 t 7 tg tg S7 S 6 S 5 Si S3 S 2 Si 

could give a nonvanishing contribution to the Dyson expansion. To characterize 
such diagrams, we begin as before by specifying in increasing time order the num- 
bers n, . . . ,r p of vertices connected through contractions within the s-block, and 
specifying the numbers h,. . . ,l q of vertices connected through contractions within 
the f-block, also in increasing time order. For example, the nonvanishing diagram 
above is described by the sequences r = (1,4,2) and I = (3,1,1,1,1,2). This 
specifies completely the (time-consecutive) contractions within the s- and i-blocks. 

It remains to specify the contractions between s- and i-variables. Note that 
we can only get additional contractions between the left endpoint of a connected 
component in the s-block with the right endpoint of a connected component in the 
t-block. Let us write m = 1 if the (left endpoint of the) zth connected component 
in the s-block is contracted with a vertex in the t-block, and «j = otherwise; 
similarly, we write Aj = 1 if the (right endpoint of the) ith connected component in 
the i-block is contracted with a vertex in the s-block, and Aj = otherwise (note 
that necessarily = S-M- For example, the nonvanishing diagram above is 
described by n — (0,1,1) and A = (1,0,0,1,0,0). Finally, we denote by n(i) the 
ith nonzero element of k, and similarly for A(i). For example, in the nonvanishing 
diagram above, k(1) — 2, k(2) — 3, and A(l) = 1, A(2) = 4. Once we have 
given r, I, k and A we have described uniquely one nonvanishing diagram, as the 
order in which the s-t contractions are made is fixed by the requirement that the 
corresponding lines be noncrossing (we must connect the lines from the inside out, 
i.e. connected component is contracted with connected component X(i)). 

With this somewhat tedious notation, we can write out the limiting Dyson series 
explicitly. Applying Wick's lemma, retaining only the nonvanishing diagrams, and 
taking the limit as e — > + gives 

oo oc £ r=n £ l=m ,_^„_ m £ «=£ A 

EEE E E 2 n-n2r»-m E / dt A ■ ■ ■ dt, d Sfl ■ ■ ■ d Sl 

B=0m=0»,mri,...,r s h, — ,l* k,\ J A rh{V J 

x (V-iI^a^i) • ■■E^MXe!${8 A ) ■ ■■Ej${8 1 )ik) l[S(s K{i) - t A(i) ), 

i 

where we have written 

pW m = / E aP {t) r = 1, 

af>{) 1 E al (t)(E n (t)y- 2 E 1() (t) r>2. 

We could proceed at this point to resum the Dyson series as before, but instead it 
will be more convenient to work backwards from the desired result and show that 
we can recover the expression above. 

Consider once more the Hudson-Parthasarathy equation 

dU t - {in dA t + L w dA\ + Lqi dA t + L 00 dt} U t . 
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We are interested in the matrix element 

(il>i\U}XU t ih) = (ai|a 2 )a«c<«i ® $|e A(/l) u}xU t e A ^ v 2 ® $). 
Using the Ito rules, we can commute the field operators past the unitaries; then 

— (ai|Q!2)o8c(/l|/2)resv(ui ® 3>|t/ t + X£7 t 1>2 <g> <&), 

where we have written 

dU t = {in(t) rfA t + Lio(i) dA\ + L 01 (t) dA t + L 00 (t) dt} U u 
dU+ = U+{LiS)d^t + L+ Q {t)dA\ + L+S)dA t + L+ {t)dt}, 
and where the coefficients are given by 

L n (t) = L n , L 10 (t)=L w + L n f2(t), L 01 (t) = L i + fi{t)*Lu, 
Lti(t) = Ll, Lt (t)=Ll 1+ L\j 2 (t), L+ l (t) = L\ + h(tYL\ 1 , 

and 

Loo(t) = £[/i(t)TW/2Wf, L+ (t) = ^[/i(t)1 a Ltj/ 2 (t)]/3. 

a/3 a/3 

But by explicit summation one may verify that 

f (t\- f + m _ y gS^) 



r>l 



-* (-i) r 2 r 
>i y ' 



Using Picard iteration to develop Ut and C/ t + into their chaotic expansions, substi- 
tuting the above expressions for L, L + and rearranging the summations somewhat, 
we arrive at the following Dyson expansion for (vi ® <fr\U^~ XU t v 2 <8> 

oo oo E r =™ E '="i , -\n-m 

EEE E E E-E E-E 

n=0 m=0 ft,m ri,...,rft ii, •••,(«! anPfi ai/3i Mm"* 



xXx f &™ (s A ) ■ ■ ■ &<£l ( Sl ) dA«f» ■ ■ ■ dA^ v 2 ® $) . 



Using the quantum Ito rules and by induction on the iterated integrals, it is not 
difficult to establish that the vacuum matrix element in this expression vanishes if 
any of the /ii or Pi are nonzero, or if the number of nonzero i/s and a's do not 
coincide. Hence we find, relabeling the variables suggestively, 

oo oo E r = n E l=rn . - n _ m E «=E A 

yyy y y — ^ y 

n—0 m—0 n,m n,...,rft k,A 



xXx / ^ >( aft )...^( ai )dA5f°...dA^%®$). 
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But now we can easily reduce to the previous form of the Dyson expansion, taking 
into account the identity (which follows directly from the quantum Ito rules) 

/ F T dA T x I G a dA ] a w®<&) 
Jo Jo 

= I dr f d<j(v(g>$\F T G a w(g)$) 5(t - a). 
Jo Jo 

The proof of Theorem 2 is complete. 

6.3. Proof of Theorem 15.11 The hard work has already been done in the proof 
of Theorem 14.31 all we have to do to prove Theorem 15.11 is an appropriate shift of 
the coefficients. We briefly provide the details. Consider first the expansion for 
(^\U t (e^W(g t] )U t (e)^ 2 ), 



oc oc „ „ 

EEH) n ~ m / dt m ---dt x 

n=0 m =0 -/A m (t) J A 



'A„(t) 



x W(g t] - 2. 9 + JB( f g+ (0, E ))T,„ (s) ■ ■ ■ T S1 (e) v 2 ® $) , 

where we have dropped the prefactor (ai|o!2)osc (/i|/2)resv and the constant factor 
that is obtained from commuting e A<y ^ 1 \ etc., past the Weyl operators. Splitting 
up the Weyl operators as explained in section [3] and commuting them through the 
Hamiltonians T as in lemma PTTI gives 



dt r , 



■dtt 



ri—0 m—0 



A m (t) 



ds n ■ ■ ■ dsi 

/A„(t) 

x ( Vl ® • • • T^( £ )f 1(e) ■ ■ ■ tl(e)v 2 ® $>, 

where we have dropped the constant factor that is obtained when we split the Weyl 
operators. Here is obtained from T s by transforming 

4e 

{st]- 2 4", e } o> £ )- 



J5io(s,e) >-> J5io(s,e) 
£oo(s,e) >-> ^oo(s,e) 



ff +(0, e) G £ (s) 



[9t] ~ 2 3 +J-( s , e) - - 3 +(0, e) G £ (s) 



and is obtained from T s by transforming 



E Q1 (s,e) h-> i?oi(s,e) 



^oo(s,e) i-> £bo(s,e) 



As 



{ gt] -2gl e }-(s,ey--g+(0,eyG e (s) 



1 
4s 



to] -24 e }-( S , £ )*--4(0, E )*G E ( S ) 



£n(s,e), 
£m(s,e), 

£n(s,e), 
^io(s,e). 



We can now proceed exactly as in the proof of Theorem 14. 31 to establish that in the 
limit e — > + , this expansion reduces to 



J2 T—n l—fn 



EEE E E 

n— m—0 n,m ri,...,r^ d,...,^ 



(-»)»- 



>n— nom- m 



E 



dtjYi • • • dt\ 



A A (t) 

K r i) v / 



dsft • • • dsi 



A fl (i) 



re(i) 



*A(i)) 5 
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where E^l (s) is obtained through the replacements 

E w (s) h-> E 1Q (s) - g(s)E 11 (s), E 00 (s) h-> E 00 (s) - g(s)E 01 (s), 
and E^q (s) is obtained through the replacements 

Eoi(s) i * Boi(s) + s(s)*£n(s), ^oo(s) i ^ Soo(s) +5(s)*^io(s). 
Starting from the opposite direction, it is not difficult to establish that the expected 
result of Theorem 15. II (tpi\U t ;W(— gt^Ut^), can be written (modulo prefactor) as 
(vi ® $|t/ t A [/ 4 v V2 ® where !7 t A is obtained from C7 t + by the replacements 

L+( S ) £+(«) - fl ( a )i+(s), Z+ (*) -> £+ (,s) - .g( S )L+( S ), 

and f/ t v is obtained from [/* by the replacements 

Lqi(s) h-> Lqi(s) +g(s)*Ln(s), £oo(s) h-> Loo(s) + g(s)*L w (s). 

It is important to note that the constant factor which we have dropped here is 
precisely the limit as e — ► + of the constant factor that was dropped previously; 
hence it suffices to show that the two expansions above coincide. However, this is 
immediate from our previous results, and the theorem is proved. 

Acknowledgments. The authors thank Luc Bouten and Hideo Mabuchi for insight- 
ful discussions. R.v.H. thanks Andrew Doherty and Howard Wiseman for their 
comments on the relation of this paper to previous work. 

References 

[1] Accardi, L., Frigerio, A., Lu, Y.: The weak coupling limit as a quantum functional central 

limit. Commun. Math. Phys. 131, 537-570 (1990) 
[2] Accardi, L., Gough, J., Lu, Y.: On the stochastic limit for quantum theory. Rep. Math. Phys. 

36, 155-187 (1995) 

[3] Barchielli, A.: Continual measurements in quantum mechanics and quantum stochastic calcu- 
lus. In: S. Attal, A. Joye, C.A. Pillet (eds.) Open Quantum Systems III: Recent Developments, 
pp. 207-292. Springer (2006) 

[4] Belavkin, V.P.: Quantum stochastic calculus and quantum nonlinear filtering. J. Multivar. 
Anal. 42, 171-201 (1992) 

[5] Billingsley, P.: Convergence of probability measures. John Wiley & Sons Inc., New York 
(1968) 

[6] Bouten, L.M., Van Handel, R., James, M.R.: An introduction to quantum filtering, preprint, 

|arXiv:math.OC/060174T1 (2006) 
[7] Doherty, A.C., Parkins, A.S., Tan, S.M., Walls, D.F.: Motional states of atoms in cavity 

QED. Phys. Rev. A 57, 4804-4817 (1998) 
[8] Dunningham, J. A., Wiseman, H.M., Walls, D.F.: Manipulating the motion of a single atom 

in a standing wave via feedback. Phys. Rev. A 55, 1398-1411 (1997) 
[9] Fagnola, F.: On quantum stochastic differential equations with unbounded coefficients. 

Probab. Th. Rel. Fields 86, 501-516 (1990) 
[10] Gardiner, C.W.: Adiabatic elimination in stochastic systems. I. Formulation of methods and 

application to few-variable systems. Phys. Rev. A 29, 2814-2822 (1984) 
[11] Gardiner, C.W., Zoller, P.: Quantum Noise. Springer, Berlin (2000) 

[12] Gough, J.: A new approach to non-commutative white noise analysis. C. R. Acad. Sci. Paris 
Ser. I Math. 326, 981-985 (1998) 

[13] Gough, J.: Quantum white noises and the master equation for Gaussian reference states. 
Russ. J. Math. Phys. 10, 142-148 (2003) 

[14] Gough, J.: Quantum flows as Markovian limit of emission, absorption and scattering inter- 
actions. Commun. Math. Phys. 254, 489-512 (2005) 

[15] Gough, J.: Quantum Stratonovich stochastic calculus and the quantum Wong-Zakai theorem, 
preprint, arXiv:math-ph/0511046 (2006) 



SINGULAR. PERTURBATIONS OF QSDE 



23 



[16] Hudson, R.L., Parthasarathy, K.R.: Quantum Ito's formula and stochastic evolutions. Com- 

mun. Math. Phys. 93, 301-323 (1984) 
[17] Marcus, S.: Modeling and analysis of stochastic differential equations driven by point pro- 
cesses. IEEE Trans. Inf. Th. 24, 164-172 (1978) 
[18] McShane, E.J.: Stochastic calculus and stochastic models. Academic Press (1974) 
[19] Meyer, P.A.: Quantum Probability for Probabilists. Springer, Berlin (1993) 
[20] Parthasarathy, K.R.: An Introduction to Quantum Stochastic Calculus. Birkhauser, Basel 
(1992) 

[21] Schweber, S.S.: An introduction to relativistic quantum field theory. Dover, Mineola, New 
York (2005) 

[22] Thomson, L.K., Mancini, S., Wiseman, H.M.: Continuous quantum nondcmolition feedback 

and unconditional atomic spin squeezing. J. Phys. B: At. Mol. Opt. Phys. 35, 4937 (2002) 
[23] Verhulst, F.: Methods And Applications Of Singular Perturbations. Springer (2005) 
[24] Warszawski, P., Wiseman, H.M.: Adiabatic elimination in compound quantum systems with 

feedback. Phys. Rev. A 63, 013,803 (2001) 
[25] Wiseman, H.M.: Quantum theory of continuous feedback. Phys. Rev. A 49, 2133-2150 (1994) 
[26] Wiseman, H.M., Milburn, G.J.: Quantum theory of field-quadrature measurements. 

Phys. Rev. A 47, 642-662 (1993) 
[27] Wong, K.S., Collett, M.J. , Walls, D.F.: Atomic juggling using feedback. Opt. Commun. 137, 

269-275 (1997) 

John Cough, School of Computing and Informatics, Nottingham Trent University, 
Nottingham, NCI 4BU, UK 

E-mail address: j ohn . goughOntu . ac . uk 

Ramon van Handel, Physical Measurement and Control 266-33, California Institute 
of Technology, Pasadena, CA 91125, USA 
E-mail address: ramonaits.caltech.edu 



